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FRONTISPIECE 


The  prototype  buoy  system  which  stimulated  this  project  is 
shown  above.  This  buoy,  part  of  the  Hysurch  Program  of  the  U.  S. 
Naval  Oceanographic  Office,  served  as  the  model  for  the  analysis 
reported  herein. 


A COMPUTER  STUDY  OF  THE  HYDRODYNAMIC  LOADINGS 
OF  A COMPLIANT  TAUT-WIRE  MOORING 

ABSTRACT 


\ 


This  report  presents  the  results  of  a computer  study  to  inves- 
tigate, primarily,  the  effects  of  hydrodynamic  loadings  on  a highly 
compliant  taut-vire  buoy  mooring  system  A computer  simulation  of  a 
prototype  buoy  system  was  used  to  study  the  influence  of  ocean  currents 
(0  to  2.8  knots),  wind  loads  (0  to  30  knots),  water  depth  (down  to  ^86  ft.) 
and  mooring  compliance  on  the  behavior  of  a highly  elastic  mooring. 
Prototype  experiments  on  the  buoy  studied  herein  have  been  conducted 
and  reported  elsewhere.^ 
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superscript  notation  is  used  to  denote  references  listed  in  Appendix  I 


v 


TABLE  OF  CONTENTS 


ACKNOWLEDGEMENTS 

ABSTRACT 

LIST  OF  ILLUSTRATIONS 

I.  INTRODUCTION 

A.  Background 

B.  Objectives  of  this  Report 

II.  BACKGROUND  AND  ASSUMPTIONS  USED  IN  COMPUTER  MODEL 

A.  Cable 

1.  Hydrodynamic  Drag  on  Cable 

2.  Cable  Weight  and  Buoyancy  Forces 

B . Surface  Buoy 

1.  Buoy  Volume  and  Projected  Area  Normal  to  Current 

2.  Hydrodynamic  Forces  on  the  Buoy 

3.  Buoy  Weight  and  Buoyancy  forces 

C.  Parabolic  Weights 

D.  The  Tripod 

1.  Tripod  Drag  Due  to  Current 

2.  Weight  and  Buoyancy  Forces  of  the  Tripod 

E.  Antenna  and  Base 

1.  Forces  Due  to  Wind  Loading 

2.  Weight  Force  of  the  Antenna 

F.  Computer  Program 

1.  General  Procedure 

2.  Node  and  Segment  Length  Selection 

3.  Velocity  Profile  of  Ocean  Current 

t 


Page 

iii 

v 

ix 

1 

1 

k 

5 

5 

8 

10 

10 

10 

11 

11 

12 

12 

12 

14 

14 

Ik 

15 
15 

15 

16 
18 


f 

•* 

5 


' 


Table  of  Contests  (cont.) 


Page 


III.  RESULTS  OF  COMPUTER  SIMULATIONS 

A.  Current  Versus  Horizontal  Displacement 

B.  Effect  of  Wind  Loads  on  Horizontal  Displacement 

C.  Link  Tension  As  a Function  of  Horizontal  Displacement 

D.  Horizontal  Displacement  as  a Function  of  Water  Depth 

E.  Compliance  of  Parallel  Links 

F.  Comparison  of  Simulation  With  Experiment 
IV.  SUMMARY  AND  RECOMMENDATIONS 

APPENDIX  I REFERENCES 

APPENDIX  II  DESIGN  DATA  FOR  PROTOTYPE  BUOY  USED  IN  THE 
COMPUTER  STUDY 

APPENDIX  III  COMPUTER  PROGRAM  LISTING 


LIST  OF  ILLUSTRATIONS 


FIGURE 


Lumped  Parameter  Representation  of  a Cable 


3 Configuration  of  the  Tripod 


Crosa-Sectional  View  of  Surface  Buoy  and  Anchor 


With  Nodal  Point  Comparisons 


5 Compliance  of  One-Inch  Rubber  Links 


Displacement  of  a Compliant  Taut-Wire  Mooring 
(0-1.6  Knots) 


Displacement  of  a Compliant  Taut-Wire  Mooring 


(1.6-2. 8 Knots) 


Displacement  for  a Compliant  Taut-Wire  Mooring 


Displacement  of  a Compliant  Taut-Wire  Mooring 


Effect  of  Wind  Load  Alone  on  the  Surface  Buoy 


Displacement  for  a Compliant  Taut-Wire  Mooring 


Effect  of  Ocean  Current  on  Mooring  Line  Tension 


in  a Compliant  Taut-Wire  Mooring 


12  Effect  of  Mooring  Depth  on  Horizontal  Displacement 
of  Compliant  Taut-Wire  Mooring  (0.8  Knot  Current) 

13  Effect  of  Mooring  Depth  on  Horizontal  Displacement 
of  a Compliant  Taut-Wire  Mooring  (1.6  Knot  Current) 


i ir  it  i ni— mi 


List  of  Illustrations  (coat.) 

Figure 

Page 

Ik 

Effect  of  Mooring  Depth  on  Surface  Buoy  Displacements 

For  a Compliant  Taut-Wire  Mooring  for  Two  Ocean  Current 

Velocities 

31 

15 

Effect  of  Ocean  Currents  on  the  Horlsontal 

Displacements  of  a Compliant  Taut-Wire  Mooring 

32 

16 

Effect  of  Mooring  Line  Compliance  on  Horisontal 

Displacement  for  Three  Different  Compliant  Moorings 

3k 

APPENDIX 

I 

Mooring  Details  for  Shallow  Water  ( 52 ' ) Installation 

42 

IX 

Mooring  Details  for  Deeper  Water  (lUO • ) Installation 

46 

I.  INTRODUCTION 


A.  BACKGROUND 

In  February  1967,  the  Experimental  Astronomy  Laboratory  of 
the  Massachusetts  Institute  of  Technology  commenced  a systems  analysis 
and  feasibility  study  of  a new  generation  hydrographic  survey  and 
charting  system  for  the  U.  S.  Naval  Oceanographic  Office.  As  a con- 
sequence of  that  study,  published  in  February  19681,  a new  approach 
for  hydrographic  surveys  of  coastal  regions  was  proposed.  The  proposed 
system,  known  as  HYSURCH  (Hydrographic  Survey  and  Charting  System) 
includes  all  aspects  of  a hydrographic  survey,  ranging  from  high  speed 
navigation  and  bathymetric  data  collection  to  rapid  production  of 
completed  charts  of  the  surveyed  areas.  One  of  the  basic  purposes  of 
the  HYSURCH  System  is  to  decrease  the  time  required  to  survey  a new 
area  and  to  produce  and  publish  charts.  (For  a detailed  description 
of  this  new  system,  the  reader  is  referred  to  Reference  1.) 

The  HYSURCH  System  uses  radio-navigation  buoys,  of  the  range- 
range  or  hyperbolic  type,  for  the  navigation  network.  Bathymetric 
data  is  acquired  aboard  high  speed  (up  to  UO  knots)  survey  vessels. 

An  error  analysis  ^ of  positional  accuracy  of  bathymetric  data  indi- 
cates that  navigation  buoy  locations  must  be  known  within  specific 
limits.  Hence,  the  "watch"  circle  specifications  for  the  navigation 
buoys  have  been  set.  To  determine  the  performance  characteristics  of 


navigation  buoys,  an  engineering  design  study  was  conducted  by  the 


Ubiversity  of  Nev  Hampshire.  Various  buoy  configurations  vere  ana- 


lysed, and  a particular  design  approach  vas  chosen  for  full-scale 


prototype  ev*'  at ion  tests.  The  results  of  this  study,  conducted  in 


various  ocean  current  and  sea  states,  and  in  several  vater  depths  are 


published 


Concurrent  vith  the  field  studies,  analytical  models  vere 


developed  to  predict  "watch"  circle  performance  as  a function  of  ocean 


currents,  vind,  vater  depth,  and  basic  design  parameters.  In  July 

, U 

1969,  an  internal  report  vas  submitted  to  the  University  of  Nev 


Hampshire  in  vhich  a computer  analysis  of  the  buoy  system  vas  conducted 


and  the  results  compared  to  field  measurements . This  computer  program 


written  by  Pladars,  is  based  upon  the  work  of  Skop,  O'Hara  and  Kaplan 


The  buoy  configuration  studied  in  the  full-scale  prototype 


experience^  and  in  the  computer  model  of  Pladars,  utilizes  a stabilized 
parabolic  surface  buoy  (see  frontispiece)  and  highly  compliant  or 


elastic  mooring  system.  A detailed  description  of  this  approach  is 


contained  in  Reference  3 


The  elastic  mooring  technique  is  not  a nev  approach  to  buoy 


technology.  The  elastic  characteristics  of  non-metallic  mooring  lines 


have  long  been  considered  important  to  mooring  performance.  One  of 
the  early  uses  of  elastic  mooring  lines  (actually  rubber  stock  cord) 


vas  in  the  surface  marker  buoys  in  the  Mohole  Project  - Phase  I,  in 


The  rubber  strands  have  been  tested  in  a buoy  installation  in 


Los  Angeles  harbor  by  D.  E.  Maddux  of  Marine  Advisors  Inc.,  during 


The  rubber  links  in  the  Marine  Advisors  Inc.  installation  evidenced  little 

» 

deterioration  during  lU  months  of  exposure  to  sea  water.  G.  H.  Savage 
and  J.  B.  Hersey,  in  the  umbilical  cord  design  for  the  Seaspider  Buoy 

O 

System,  used  four  1"  diameter  35  ft.  solid  rubber  strands  that  form  a 
compliant  connection  between  surface  and  subsurface  buoys.  John  Isaacs, 
of  the  Scripps  Institution  of  Oceanography,  has  employed  compliant 
mooring  to  deep-ocean  systems The  basic  principle  in  all 
these  moorings  is  the  use  of  elastic  materials  with  high  compliance 
to  reduce  the  changes  in  mooring  tensions  due  to  relative  displacements 
in  the  moorings 'lines , i.  e.  waves,  tides,  and  lateral  movement  by  currents. 

Out  of  the  in  situ  experience  of  these  compliant  mooring  designs, 
it  became  increasingly  obvious  that  some  mathematical  modelling  of 
the  approach  should  be  conducted.  As  an  initial  step  in  modelling 
such  systems,  Pladars**  studied  the  performance  of  a specific  iesign 
configuration  to  steady  state  ocean  currents  and  winds.  The  effect 
of  wave  action  per  se  was  neglected,  because  the  particular  buoy  system 
specification1  included  substantial  ocean  currents  but  only  modest 
sea  conditions  (up  through  Sea  State  U).  As  a result  of  Pladars'  study, 
a computer  model  existed.  Therefore,  it  was  deemed  desirable  to 
continue  to  utilize  it  to  study  the  effects  of  changing  certain 
environmental  conditions  and  design  parameters.  This  report  is  con- 
cerned with  reporting  the  results  of  this  latter  study. 


•No  reference  appears  to  be  available.  Statements  reported 
herein  came  from  informal  communications. 
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This  report  presents  the  result  of  a computer  study  of  the 
performance  characteristics  of  the  HYSURCH  Navigation  Buoy  discussed 
in  Reference  3.  The  computer  program  (listed  in  Appendix  III), 
written  by  Pladars^  in  Fortran  IV  for  an  IBM  36O/UO  system,  is 
used  to  study  the  relation  between  steady  state  displacements  and 
configurations  when  variables  such  as  ocean  current,  wind  velocity, 
depth,  rubber  link  tension  and  compliance  were  altered  from  a given 
design^  configuration  of  a taut-wire  surface  moored  parabolic  buoy. 

The  results  reported  herein  are  limited  to  a parabolic  buoy  system 
employing  an  elastic  link  as  part  of  the  mooring  design.  The  problem 
is  considered  to  be  two-dimensional.  The  study  should  be  considered  as 
preliminary,  for  many  more  variables  still  need  to  be  considered. 

It  does  , however,  present  some  trends  in  the  behavior  of  compliant 
moorings . 

The  report  is  divided  in  parts , as  follows: 

Section  I - Introduction 

Section  II  - A brief  discussion  of  the  assumptions 
used  in  the  computer  program 

Section  III  - A discussion  of  the  results  of  the 
computer  simulation. 

Section  IV  - Summary  and  Recommendations 
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II.  BACKGROUND  AND  ASSUMPTIONS  USED  IN  COMPUTER  MODEL 


I 


A relatively  simple  computer  model  was  developed,  based  on 
these  types  of  assumptions: 

1)  Hydrodynamic  forces  are  modelled  using  lumped 
parameter  techniques. 

2)  Wind  induced  forces  on  the  buoy  structures  above 
water  are  simulated  using  the  "fifth  power  law" 
for  the  velocity  profile  and  by  lumped  parameter 
representation . 

3)  All  aero /hydrodynamic  forces  acting  on  the  buoy 
system  are  viewed  as  coplanar. 

U)  Wave  generated  forces  are  neglected. 

This  section  of  the  report  describes  the  various  aspects  of  the 
computer  model,  including  a summary  of  the  mathematical  assumptions 
made.  The  specific  design  parameters  are  listed  in  Appendix  II,  with 
a schematic  diagram  shown  in  Figure  1. 

A.  CABLE 

For  the  cable  analysis,  a lumped  parameter  representation  of 

1 P 

the  cable  was  employed,  similar  to  that  used  by  Paquette  and  Henderson. 

In  a lumped  parameter  representation  the  cable  is  modelled  as  a series 
of  stations  or  nodes,  connected  by  straight-line  segments  of  cable 
(see  Figure  2).  All  forces  due  to  current,  weight,  and  buoyancy  of 


' 
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the  cable  half-segment  adjacent  above  and  below  the  node  are  assumed 
to  be  concentrated  at  each  particular  node.  Ths  cable  segments  func 
tlon  only  to  connect  the  nodes  and  are  assumed  to  present  no  resis- 
tance to  bending. 

1.  Hydrodynamic  Drag  on  Cable 

/ 

Only  drag  forces  normal  to  the  cable  verb  considered.  It  was 
felt  that  tangential  drag  vould  be  of  little  significance  for  this 
particular  system  since  the  cable  lengths  are  less  than  500  feet  and 
because  the  orientation  is  in  most  cases  virtually  vertical. 

Using  the  cross-flow  principle  for  inclined  circular  cylinders 
and  the  fluid-dynamic  drag  relationship, ^ the  current  drag  force 
normal  to  the  cable  segment  can  be  written, 

PU  - 1/2  p CD  AV2  cos2  0 Eq.  (l) 

Fn  ■ 1/2  p CD  dlV2  cos2  0 Eq.  (2) 

yhere:  Fjj  ■ Hydrodynamic  Force  Normal  to  Cable 
p ■ Density  of  Water 
CD  = Drag  Coefficient 

A = Area  of  Cable  Normal  to  Local  Fluid  Velocity 
d = Cable  Diameter 

1 = Length  of  Cable  Over  Which  FN  Acts 

V ■ Local  Fluid  Velocity 

6 * Angle  between  Cable  and  Local  Fluid  Velocity 


2 . Cable  Weight  and  Buoyancy  Forces 


Due  to  the  constant  volume  assuaption  for  the  rubber  links, 
the  cable  weight  and  buoyancy  forces  can  be  computed  on  the  basis  of 
the  original  unstressed  dimensions: 

PV  " WCAB  10  Eq.  (5) 

Ffi  • pgwd2l/l»  Eq.  (6) 

B.  SURFACE  BUOY 

1.  Buoy  Volume  and  Projected  Area  Normal  to  Current 

An  8-foot  parabolic  buoy  was  used  in  the  prototype  studies. 

It  is  constructed  essentially  of  two  paraboloids  of  revolution. 

Since  the  actual  paraboloid  formula  was  not  available  from  the  manufacturer 
2 

r = 12h  was  determined  to  provide  an  approximation  of  the  projected 
area  normal  to  the  current,  where  r (major  radius)  and  h (minor  radius) 
ar®  io  feet.  Thus,  the  volume  and  cross-sectional  area  of  one 
parabolic  half  can  be  computed  as: 

Vol  * 1/8  irh  (2r)2  Eq.  (7) 

= 1/8  irh  (48h) 

= 6 irh2 

Area  - 4/3  hr  Eq.  (8) 

= 8/3  /3  h 3/2 
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c Forces  on 


Assuming  that  the  drag  on  the  buoy  is  proportional  to  the 
submerged  cross-sectional  area  of  the  buoy  normal  to  the  current  and 
neglecting  any  wave  forces  and  small  inclinations  of  the  buoy  (mea- 
sured in  prototype  studies  to  be  less  than  8°),  the  drag  force  could 
be  written  as 


Since  the  tensions  of  the  cable  are  such  that  more  than  half 


area  is  simply  the  total  cross-sectional  area  minus  the  cross-sectional 


The  lift  force  on  the  buoy  was  disregarded  because  the  lift 


was  estimated  to  be  insignificant  compared  to  buoyancy  force  changes 


due  to  high  buoy  volume  concentration  at  the  waterline  for  this 


3.  Buoy  Weight  and  Buoyancy  Forces 


Weight  of  the  buoy  was  not  determined  separately  but  incorpor 


ated  into  the  total  weight  of  the  system 


Buoyancy  force  on  the  buoy  is  given  by  the  following  expression 


Buoyancy  Force 


p » Mass  Density  of  Water 


Gravitational  Constant 


Vq  ■ Submerged  Volune 


where  the  submerged  volume  is  obtained  by  subtracting  the  volume  above 
the  waterline  from  the  total  voluae  of  the  buoy. 

C.  PARABOLIC  WEIQHTS 

The  parabolic  daaper  weights  (aee  disks  - shaped  eleaents 
at  base  of  buoy  in  Figure  l)  attached  at  the  base  of  the  bridle 
or  tripod  are  shaped  similarly  to  the  buoy.  The  determination  of 
forces  due  to  these  weights  follows  the  same  procedure  as  for  the 
buoy.  The  only  difference  is  that  the  parabolic  weights  are  fully 
submerged  at  all  times. 

D.  THE  TRIPOD 

1.  Tripod  Dreg  due  to  Current 

Vertical  forces  due  to  current  have  been  disregarded  because 
they  are  small  compared  to  other  forces  at  the  buoy. 

Leg  1 of  the  tripod  (see  Figure  3)  can  be  handled  in  exactly 
the  sane  manner  as  a cable  segment,  thus  giving 

Fx  - 1/2  p Cpdl  V2  cos3  - Eq.  (ll) 

vhere:  Fx  ■ Hydrodynamic  Force  on  the  Tripod 

* ■ Angle  between  Local  Fluid  Velocity  and  the  Tripod  Legs 
Note:  other  terms  similarly  defined  as  in  Equations  142 

For  legs  2 and  3,  from  symmetry  basic  hydrodynamic  stability  of 
the  orientation,  and  forces  perpendicular  to  the  x-y  phase  will 
cancel.  Taking,  say,  leg  2 and  resolving  the  current  into  vector 
components,  • |V|  sin  30°  and  |$g|  » |V|  sin  30°,  each 
component  can  be  handled  separately  as  an  effective  current.  The 

J 
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Fx  - 1/2  p CD<11  V*  cos  30°  Eq.  (li 

in  the  direction  of  V^,  or 

- 1/2  p Cpdl  V2  cos3  30°  Eq.  (l 

in  the  x-direction.  The  dreg  force  on  leg  2 due  to  V2,  using  the 


cross-flov  principle,  becomes 


in  the  direction  of  V, 


F2x  ■ 1/2  p CDdl  V*  sin-1  30°  cosJ  « Eq.  (15) 

in  the  x-direction.  Combining  the  x-direction  forces  for  all  three 


Veight  of  the  tripod  is  included  in  the  weight  of  the  system 


Buoyancy  is,  again,  simply  a function  of  the  volume,  yielding  an 


expression  similar  to  a cable  segment 


1.  Forces  due  to  Wind  Loading 


For  simplicity  the  antenna  was  assumed  to  be  vertical  at  all 

times  (measured  in  prototype  studies  to  be  less  than  8°  off  vertical 
for  all  conditions  reported3).  The  antenna  vas  considered  to 


be  a vertical  rod  without  any  guy  vires  or  other  protuberances.  The 
wind  was  assumed  to  vary  with  height,  according  to  the  Fifth  Power 
Lav,  as  given  by  Saunders ,^5  although  the  logarithmic  distribution 


could  have  equally  as  well  been  used. 

\ * v6  <|>1/5 


Eq.  (IT) 


where:  is  the  wind  velocity  at  any  height  (h) 

Vg  is  the  velocity  of  wind  observed  at  6 feet 
Using  the  fluid-dynamic  drag  relationship, 

Fc  - 1/2  p CD  A V2  Eq.  (18) 

and  p ■ 0.00238,  C = 1.2,  an  empirical  formula  for  antenna  drag  was 
calculated  to  be 

Fm  * 0.011221  V2  Eq.  (19) 

and  for  the  antenna  base 

Fdb  - 0.00171  V2  Eq.  (20) 

giving  a combined  total  drag  of 

Fd  » 0.012931  V2  Eq.  (21) 

2 . Weight  Force  of  the  Antenna 

Weight  of  the  antenna  and  attached  equipment  was  included  in 
the  total  weight  of  the  system. 

F.  COMPUTER  PROGRAM 
1.  General  Procedure 

Initially,  all  forces  acting  on  the  antenna,  the  buoy  the 
tripod,  the  parabolic  weights,  and  the  top  half-segment  of  the  cable 


•re  euaned  up  end  applied  to  the  top  end  of  the  cable.  Starting  from 
this  top  end  of  the  cable  and  proceeding  tovards  the  anchor  end,  the 
resultant  forces,  tensions,  new  dimensions , and  orientations  of  each 
cable  segment  are  determined.  This,  in  turn,  locates  the  nev  position 
of  each  station  or  node.  When  the  boundary  conditions  for  the  top  end 
of  the  cable  are  satisfied,  all  forces  acting  on  the  top  and  along  the 
cable  are  recomputed.  This  iteration  process  continues  until  the  node 
coordinates  between  successive  computations  change  less  than  the  speci- 
fied value.  The  computation  is  then  considered  complete  and  desired 
output  is  printed. 

2.  Node  and  Segment  Length  Selection 

In  the  selection  of  segment  lengths  some  consideration  must  be 
given  to  the  expected  radius  of  curvature  of  the  cable:  a small  radius 
of  curvature  should  have  nodes  closer  together  them  a large  radius  of 
curvature.  A node  should  also  be  located  at  every  point  where  a prop- 
erty of  the  cable  changes  or  an  external  element  is  present.  A study 
of  the  effects  of  numbers  of  nodes  was  made.  The  number  of  nodes  se- 
lected is  a function  of  the  particular  situation  being  studied,  although 
15  nodes  were  used  in  all  the  calculations.  The  effect  of  increasing 
the  number  of  nodes  from  15  to  22,  for  example,  is  noted  on  Figure  H 
(vhich  indicates  that  the  maximum  difference  in  the  horizontal  displace 
ment  of  any  nodal  point  is  O.OlX  of  these  two  cases).  Other  studies 
were  conducted  to  determine  effects  of  node  selection. 


Figure  4 Croaa  Sectional  View  of  Surface  Buoy  and  Anchor 
with  Nodal  point  Coawarlaona 


3.  Velocity  Profile  of  Ocean  Current 


i. 


Since  the  design  of  the  Bystem  vas  originally  for  shallov 
water,  a constant  surface  to  bottom  profile  vas  selected.  Use  of 
the  program  for  other  than  constant  current  profile  will  require  . 

minor  modification  of  the  program.  Skop  and  Kaplan, ^ for  example, 
consider  a velocity  profile  that  can  be  represented  by  straight-line 
segments . 

III.  Results  of  Computer  Simulations 

This  section  of  the  report  contains  the  results  of  several 
computer  simulations  to  determine  the  effects  on  system  behavior 
due  to  changes  in  environmental  conditions  or  changes  in  design 
parameters . 

4 

A.  CURRENT  VERSUS  HORIZONTAL  DISPLACEMENT 

The  first  phase  of  the  analysis  was  to  predict  the  relationship 
between  ocean  current  and  the  horizontal,  displacement  of  a given  buoy 
configuration.  The  anchor  was  placed  136  feet  -below  the  water  suS-s* 
face  (since  experimental  data  is  available  for  this  depth).  The 
compliance  of  the  four  parallel  one-inch  rubber  links  was  expressed 
as  Case  3 in  Figure  5.  Ocean  currents  were  varied  from  0 to  2.8  knots 
at  0.4  knot  increments.  A wind  of  10  knots  was  assumed  in  the  direction 
of  the  current.  Displacement  was  plotted  as  a function  of  depth  for 
various  currents  giving  the  shape  of  the  cable.  The  results  can  be 
seen  graphically  in  Figures  6 and  7.  . 

The  shape  of  the  mooring  line  is  almost  linear  at  very  small 
currents ; hovever,  as  the  current  increases  the  mooring  line  develops 
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iom  curvature.  The  effects  of  current  on  the  displacement  of  the  surface 
buoy  is  shown  in  Figure  8. 


The  vind  velocity,  as  shown  in  the  next  section,  has  a negligible 
effect  on  this  result.  The  major  contributor  to  the  displacement  is 
the  drag  force  on  the  tripod.  It  is  calculated  to  be  more  than  twice  the 
buoy  drag  force.  The  parabolic  weights  are  less  significant  in  drag 
considerations.  However,  at  small  currents,  0 to  0.4  knots,  the  drag 
forces  on  the  tripod,  buoy  and  weights  are  completely  overshadowed  by  the 
antenna  drag  forces  due  to  wind  if  there  is  a significant  wind  (10 
knots  or  more). 

B.  EFFECT  OF  WIND  LOADS  ON  HORIZONTAL  DISPLACEMENT 

The  next  simulation  considered  the  effect  of  wind  loads  on 
the  horisontal  displacement  of  the  buoy.  An  ocean  current  of  zero 
was  assumed.  The  wind  velocity  was  taken  to  be  5,  10,  17  and  30  knots 
respectively.  The  results  can  be  seen  in  Figure  9.  The  mooring  line  node 
loci  are  linear,  which  is  to  be  expected  having  only  a wind  induced  force 
at  the  top  of  the  cable  parallel  to  the  surface. 

The  relationship  between  the  change  in  wind  velocity  versus  dis- 
placement is  similar  in  form  to  that  of  the  change  in  current  versus  the 
change  in  displacement  because  the  basic  hydrodynamics  are  the  same. 

At  wind  loads  from  approximately  0 to  20  knots  the  relation  is 
parabolic,  and  at  higher  wind  velocities  it  appears  to  become  more 
linear.  Winds  up  to  40  knots  vere  studied  and  the  apparent  linear 
relation  held. 
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, It  was  postulated  that  an  Increase  in  the  horizontal  displacement 

1 

of  the  cable  vould  not  change  significantly  the  tension  in  the 


rubber  links.  Over  the  design  range  of  operation,  250  to  300  pounds, 

£ 

the  tension  change  is  predicted  to  be  small  (about  20jt  for  a 2.8 
knot  current;  10  knot  coplanar  wind)  compared  to  the  change  in 
elongation,  as  shown  in  Figure  11.  For  example,  a change  in  displacement 
of  17  feet  varies  the  tension  only  18  pounds  at  2 knots.  Up  to 
approximately  2 knots  the  tension  and  displacement  relationship  is 
approximately  one  pound  of  tension  per  one  additional  foot  of  displacement. 
Above  2 knots,  the  tension  increases  more  rapidly. 
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D.  HORIZONTAL  DISPLACEMENT  AS  A FUNCTION  OF  WATER  DEPTH 


This  buoy  vsa  designed  for  installation  in  up  to  200  feet  of 
vater.  An  investigation  of  the  effect  of  deeper  vater  depths  on 
buoy  displacement  vas  carried  out  by  adding  segments  to  the  cable 
node  in  the  program.  The  analysis  vas  done  for  tvo  ocean  currents: 

8 knots  and  1.6  knots.  The  vater  depth  vas  increased  from  136  feet 
to  U86  feet  in  50  foot  increments.  The  results  are  shovn  graphically 
in  Figures  12  and  13.  The  change  in  horizontal  displacement  for 
each  change  in  depth  is  shovn  in  Figure  lb  vhich  is  a plot  of  the 
vater  depth  versus  horizontal  displacement  of  the  surface  buoy.  At 
a current  of  0.8  knots,  the  horizontal  displacement  increases  approximately 
2.5  feet  for  every  increase  of  50  feet  in  vater  depth.  At  1.6  knots 
the  change  is  approximately  7.6  feet  for  every  50  feet,  although  the 
linear  assumption  is  not  completely  Justified. 

The  effect  of  current  at  greater  depths  vas  also  considered. 

A 10  knot  vind  vas  assumed  and  the  current  vas  varied  from  0 to  1.6 
knots  in  increments  of  O.b  knots.  We  compared  these  results  vith  the 
results  obtained  in  the  first  part  of  this  analysis  (current  versus 
displacement  at  136  feet ) . The  results  almost  exactly  paralleled  the 
results  for  shallover  installation  depth  (Figure  6).  The  current 
versus  displacement  of  cable  graph  vas  again  approximately  parabolic 
at  small  currents  and  approximately  linear  at  higher  currents 
(see  Figure  15). 
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HORIZONTAL  DISPLACEMENT  (ft) 


(WIND  LOAD  DUE  TO  K)  KNOT  WIND) 


E.  COMPLIANCE  OF  PARALLEL  LINKS 

Computer  model*  for  the  elongation  of  the  rubber  links  are  shown 
in  Figure  5 along  with  experimental  data.  The  extension  formulas 
used  are  valid  only  for  the  range  of  tensions  encountered,  250  - 300 
pounds.  The  case  3 compliance  model  was  employed  in  the  calculations 
used  throughout  this  analysis.  Curving  fitting  routines  could  be  used 
to  better  model  the  experimental  data;  however,  it  was  decided  to 
use  a linear  approximation  for  the  first  stage  of  this  analysis. 

The  program  is  currently  being  modified  to  handle  the  non-linear 
characteristic. 

In  an  effort  to  preliminarily  evaluate  the  effects  of  altered 
mooring  line  compliance,  the  displacement  characteristics  due  to 
different  compliances  for  a fixed  environmental  condition  (1.0  knot 
sea  current,  a 10  knot  coplanar  wind,  and  an  ocean  depth  of  136  feet) 
were  simulated.  Three  different  elongation  models  (Case  I,  II  & 

III  from  Figure  5)  were  placed  separately  in  the  computer  simulation. 

The  results  appear  in  the  graph  of  horizontal  displacement  for  the 
mooring  versus  depth  along  the  moor  in  the  three  different  compliance 
cases  (see  Figure  16).  The  difference  in  buoy  displacement  between 
Case  1 and  Case  3 is  significant  (approximately  150)0.  Each 
of  the  three  linear  approximations  of  the  experimental  data  represents 
an  extension  formula  for  the  range  of  tensions  actually  encountered 
in  the  prototype  experiments. 
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conducted  in  the  open  ocean  vith  sea  conditions  up  to  about 
state  3-1*  seas,  vith  measured  winds  up  to  18  knots,  surface 
currents  up  to  0.8  knots,  and  tidal  variations  of  8.3  feet 
vith  a water  depth  of  lUo  feet  (mean  low  vater).  Under  these 
conditions,  the  3o  - buoy  position  error  was  - 9.15  feet 
(calculation  based  on  all  data  recorded  for  this  site 
through  10  September  - the  12  September  data  was  not  used 
since  only  one  of  the  four  rubber  links  in  the  mooring 
vas  intact)  about  the  mean  position,  or  a 3 o - buoy 
position  vatch  circle  of  18.3  feet.  All  recorded  data  fell 
within  a vatch  circle  of  l6.h  feet.  The  antenna  vertical 
inclination  stayed  vithin  a cone  of  - 8°,  while  the  vertical 
heave  acceleration  never  exceeded  6 ft. /sec ."2 
In  general,  the  winds  were  coplanar  with  the  ocean  currents  (which 
averaged  0.58  knots  during  this  time  vith  a 3o  of  - 0.32,  or  an  upper 
3o  limit  of  0.88  knots). 


>on  of  Simulation  vith  Experiment 


In  Reference  3,  a series  of  experiments  are  reported  in  which  a 
buoy  of  a design  similar  to  that  used  in  the  computer  simulation  vas 
used.  The  results  of  the  "Open  Sea  Tests'1  are  summarised  below: 

Sea  Teats  (Isles  of  Shoals  Tests)  These  tests  were 


The  results  reported  in  Reference  3 are  statistical  in  form, 
and  hence  are  not  directly  comparable  to  the  results  from  the 
deterministic  computer  model.  However , using  the  extreme  values  of 
0.8 8 knots  current,  18  knots  of  coplaaar  wind,  and  a water  depth  of 
136  feet,  the  computer  model  predicts  a maximum  radial  displacement 
of  U.9  feet  (using  Figures  8 and  9)  whereas  the  experimental  maximum 
excursion  for  all  data  was  8.2  feet.  Thus  the  computer  simulation 
prediction  is  only  60t  of  the  measured  value.  The  comparison, 
however,  is  not  fully  Justified  since  the  in  situ  experiment  encountered 
variability  in  ocean  current  direction,  wind  waves,  and  wind,  and 
hence  the  experiment  was  conducted  on  a statistical  basis  using  a 
large  number  of  data  points  to  predict  behavior.  The  comparison 
between  computer  simulation  and  experiment  should  be  considered 
only  as  an  indication  of  the  possible  validity  of  the  model. 

Controlled  experiments  should  be  conducted  to  more  adequately 
validate  the  computer  simulation. 

Pladars\  who  reanalyzed  some  of  the  data  reported  in  Reference 
3,  found  a similar  difference  between  experimental  results  and  computer 
model. 
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IV.  SUMMARY  AND  RECOMMENDATIONS 


The  computer  simulation  reported  herein  is  not  intended  to  be 
an  exhaustive  or  general  design  study.  The  impetus  for  the  study  vas 
a specific  buoy  system  configuration,  upon  vhich  considerable 
experimentation  had  been  conducted.  The  computer  simulation  vas 
conceived  as  a way  to  predict,  in  a general  manner,  the  behavior  of 
a typical  compliant  mooring.  While  a comprehensive  study  to  validate 
the  computer  simulation  against  experiment  would  be  extremely  desirable 
the  data  obtained  during  experiment  is  not  presently  in  a form  that 
lends  itself  to  such  a study.  The  simple  comparison  vhich  is  made 
indicates  a prediction  by  computer  simulation  vhich  is  60%  of  the 
value  experimentally  determined.  A carefully  designed  experiment 
is  still  needed  to  corroborate  this  computer  model. 

Since  the  model  has  not  been  adequately  validated  against 
experiment,  the  results  contained  in  this  report  should  be  considered 
only  as  predicting  behavior  in  a general  way.  A reanalysis  of  the 
in  situ  data,  on  a point  by  point  basis,  could  possibly  lead  to  a 
more  adequate  validation  of  the  model,  but  such  is  beyond  the  scope 
of  the  vork  reported  herein. 

The  essence  of  this  report  is  contained  in  Section  III.  This 
section  discusses  the  computer  simulation  of: 

1)  Effects  of  ocean  currents  on  horizontal  displacements 
(subsection  A) 

2)  Effects  of  vind  loads  on  horizontal  displacements 
(subsection  B) 
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3)  Link  tension  as  function  of  ocean  current  (subsection  C) 

*»)  Horizontal  excursions  as  function  of  vater  depth  and 
currents  (subsection  D) 

5)  Effect  of  elastic  characteristics  on  horizontal 
displacements,  as  function  of  ocean  currents  (subsection  E) 

6)  Short  discussion  of  the  computer  simulation  as  compared 
to  experimental  data  (subsection  F) 

It  is  recommended  that  serious  thought  be  given  to  the  validation 
of  this  and  other  computer  simulations  of  buoy  moorings.  There 
are  numerous  analytical  models;  hovever,  very  limited  experimental 
corroboration  of  the  analysis.  The  prototype  experiment  in  Reference 
3 vas  designed  for  purposes  other  than  validation  of  this  computer 
model.  Hovever,  future  studies  should  include  instrumentation  which  will 
provide  data  to  correlate  theoretical  prediction  with  actual  behavior. 

It  may  be  that  laboratory  experiments  are  required  prior  to  full- 
scale  studies,  since  open  sea  tests  do  not  often  provide  controlled 


conditions. 
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APPENDIX  II 

DESIGN  DATA  FOR  PROTOTYPE  BUOY  USED  IN  THE  COMPUTER  STUDY 

The  buoy  system  selected  for  this  computer  study  is  described 
in  detail  in  Reference  3.  As  an  aid  to  understanding  the  actual  buoy 
system  used  in  this  study,  some  design  details  are  listed  below. 

The  buoy  system  details,  needed  to  understand  the  results  of 
the  computer  study,  consist  of: 

1)  Surface  Bvoy 

2)  Antenna 

3)  Mooring  Bridle 

b)  Compliant  Taut-Wire  Mooring 

5)  Anchor 
Surface  Buoy 

The  surface  buoy  consists  of  a horisontal  altitude  parabolic 
configuration,*  vhich  is  basically  eight  feet  in  diameter.  To  accommo- 
date the  instrumentation  systems,  four  12"  dia.x  32"  deep  instrument 
compartments  were  designed  into  the  buoy.  The  buoy  is  constructed  of 
aluminum  with  spun  parabolic  heads.  The  welded  buoy  is  polyurethane 
foam  filled  to  assure  buoyancy.  The  aluminum  is  treated  for  corrosion 
resistance  and  painted  with  epoxy  paint.  The  buoy,  antenna,  and  tripod 
legs  are  shown  in  Figure  I. 
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•Model  5-200  Manufactured  by  Prodelin,  Inc.,  Hightstown,  N.  J» 


Antenna 


For  this  study,  a Columbia*  Model  222  was  modified  to  produce 
a 30  foot  (effective  r.f.  length)  antenna.  This  fiberglass  antenna 
was  mounted  on  a five- foot  (5* ) aluminum  stand-off,  which  placed  the 
top  of  the  antenna  approximately  36'  above  the  ocean  surface.  The 


upper  11'  of  the  antenna  vas  not  given  additional  support;  however 


considerable  effort  was  made  to  stabilize  the  lower  parts  of  the  antenna 


Four  diamond  spreaders  were  placed  in  quadrature  between  the  antenna 


base  and  a point  18'  above  the  base  which  is  above  a known  node  on  the 


antenna.  These  spreaders  tended  to  stabilize  the  antenna  from  natural 


frequency  oscillations.  The  antenna  vas  bolted  to  the  five- foot  stand 


off,  which  in  turn  was  bolted  to  the  surface  buoy  mounting  flange. 

The  mid-point  of  the  diamond  spreader  stabilizer  (9'  from 
antenna  base)  was  supported  by  three  guy  wires  to  the  surface  buoy  (guy 


wires  were,  in  general,  3/32"  stainless  stranded  wire,  terminated  with 
r.f.  insulators  and  "nicopress"  fittings).  In  addition,  the  five  foot 
antenna  stand-off  was  supported  by  3,  3/16"  guy  wires  to  the  surface 


The  antenna  stand-off  also  supported  the  oceanographic  buoy 


Mooring  Bridle 


The  buoy  was  designed  to  operate  with  a tripod  mooring  bridle 


•Manufactured  by  Columbia  Products,  Inc.,  Columbia,  South  Carolina 
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approximately  8'2"  in  height.  The  bridle  ( see  Figure  I.  ) is  essen- 
tial to  the  operation  of  the  buoy,  since  it  provides  the  necessary 
righting  moments  for  the  specified  vertical  stability.  The  1 5/6" 
galvanised  tripod  legs  are  bolted  to  the  surface  buoy  at  one  end  and 
to  an  apex  yoke  at  the  other.  The  yoke  provides  the  mooring  point  as 
veil  as  a mounting  point  for  the  dynamic  stabilizing  weights  (shown  in 
Figure  I.  as  "Para-Weight").  These  veights  have  the  effect  of 
substantially  increasing  the  pitch/roll  moment  of  inertia,  which  is 
required  for  dynamic  stability. 

Mooring 

The  taut-vire  mooring  was  one  of  the  critical  components  in 

the  design.  Because  of  the  wide  range  of  expected  water  depths  and 

tidal  conditions,  the  mooring  must  be  capable  of  accommodating  length 

changes.  The  study  of  moorings  was  focused  on  methods  of  overcoming 

tidal  variations  while  still  maintaining  mooring  tension  tod  essentially 

a constant  buoy  vaterline.  As  a part  of  a research  program  associated 

with  the  Seaspider  effort  of  the  Office  of  Naval  Research,  the  U.  N.  H. 

design  group  had  been  studying  techniques  to  develop  quasi-constant 

tension  taut-vire  moorings.  The  first  major  application  of  the  quasi- 

constant  tension  taut-vire  mooring  was  developed  for  Seaspider  I.  The 

basic  principle  in  these  moorings  is  to  use  a spring-like  material  with 

high  compliance  compared  to  the  deflections  anticipated  in  a design. 

This  fact  can  be  easily  seen  from  the  basic  linearised  relationship  for 

a Bpring-like  material:  Force  (F)  = displacement  (x) 

compliance  (1) 
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Therefore: 

df  = dx 
F x 

which  indicates  that  the  % change  in  displacement  results  in  an  equal 
change  in  force.  Therefore,  the  length  of  the  spring-like  material 
under  load  must  be  large  compared  to  tidal  variations  (say  10  feet), 
or  the  material  must  be  highly  compliant  (note:  a constant  force 
spring  is  infinitely  compliant ) . 

There  are  numerous  ways  of  obtaining  this  compliance.  In  an 
operational  buoy  system,  one  can  imagine  a constant  tension  winch  on 
the  mooring  cable,  which  would  serve  booh  for  storing  of  the  mooring 
cables  as  well  as  providing  the  constant  tension. 

For  the  prototype  evaluation,  a relatively  simple  approach  was 
used.  The  mooring  used  four  (U)  one-inch  (l")  diameter  solid  NATSYN 
rubber  rods,  configured  as  shown  in  Figures  I & II.  For  the  load  levels 
used,  (approximately  2000  lb.)  this  rubber  produces  approximately 
2.0  lb.  per  percent  elongation  per  rubber  rod,  which  for  this  mooring  re- 
sulted in  a force  change  of  50  lb/ft  of  water  depth  change.  For  ex- 
ample, a 10-foot  tide  would  change  the  mooring  tension  by  500  In. 

This  tidal  effect  results  in  the  buoy  waterline  change  of  2 inches 
which  is  negligible. 

The  rubber  links  were  secured  to  the  anchor  with  self-aligning 
mounting  yokes  and  shackles.  A standard  underwater  swivel  was  used  to 
eliminate  torsional  rotation  problems.  All  shackles  were  of  the  seiz- 
able  type. 
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FIGURE  II 

Mooring  Details  for  Deeper  Water  (lfcO')  Installation 
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Anchoring 

The  anchors  used  were  of  the  deed-weight  type  end  were  of 
e design  similar  to  those  reported  in  Reference  8. 
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COMPUTER  PROGRAM  LISTING 

PAM  1 ~ ~ 

//  job 

LOO  DRIVE  CART  SPEC  CART  AVAIL  PHY  DRIVE 
OOOO  IDOC  1D0C  0000 

V2  MOt  ACTUAL  BK  CORF 10  IK 

//  FOR 

•LIST  SOURCE  PROGRAM 

C DISPLACEMENT  ANALYSIS  OF  A TAUT  WIRE  PARASOL IC  BUOY 

C 

C A • ANGLE  OF  TRIPOD  LEG  FROM  VERTICAL 
C ASBUOY  « SUBMERGED  CROSS-SECTIONAL  AREA  OF  BUOY 

C AT BUOY  - TOTAL  CROSS-SECTIONAL  AREA  OF  BUOY 

C 

C B ■ 30-DEGREE  ANGLE  IN  TRIPOD  GEOMETRY 
C 

C CDBUOY  - BUOY  DRAG  COEFFICIENT 

C CDCABtN)  • CABLE  SEGMENT  DRAG  COEFFICIENT  NORMAL  TO  CABLE 

C CDPRBL  " PARABOLIC  WEIGHT  ORAG  COEFFICIENT 

C CDTRPD  ■ TRIPOD  LEG  DRAG  COEFFICIENT  NORMAL  TO  LEG 

C CONST IN)  • EXPERIMENTAL  ELONGATION  CONSTANT  FOR  CABLE  SEGMENT 

C COSTHINI  - COSINE  OF  ANGLE  BETWEEN  CABLE  SEGMENT  AND  VERTICAL  AT 

C THE  N-TH  NODE 

C 

C DIN)  » DIAMETER  OF  CABLE  SEGMENT  UNDER  TENSION 
C OBUOY  • DIAMETER  OF  BUOY 

C DEIN)  ■ 12  • CD  * DIA  * LENGTH)  FOR  THE  ELEMENT 
C DELE  - CUT-OFF  VALUE  FOR  ERROR  FUNCTION 
C DELTA  ■ POSITIVE  CONVERGENCE  FACTOR 

C OELXY  - DESIRED  ACCURACY  FOR  COORDINATES 

C DEPTH  • DEPTH  OF  WATER 

C D1AIN)  • ORIGINAL  UNSTRESSED  DIAMETER  OF  CABLE  SEGMENT 
C DPRBL  - DIAMETER  OF  PARABOLIC  WEIGHTS 

C DTRPD  - DIAMETER  OF  TRIPOD  LEGS 

C 

C E • ERROR  FUNCTION 

C 

C FBBUOY  • BUOYANCY  FORCE  OF  BUOY 

C FBCAB  • BUOYANCY  FORCE  OF  TOP  HALF  SEGMENT  OF  CABLE 
C FBPRBL  • BUOYANCY  FORCE  OF  PARABOLIC  WEIGHTS 

C FBTRPD  - BUOYANCY  FORCE  OF  TRIPOO 

C FBYIN)  • BUOYANCY  FORCE  OF  CABLE  SEGMENT 

C FCYIN)  • VERTICAL  FORCE  ON  CABLE  SEGMENT  DUE  TO  CURRENT 

C FCYCAB  • VERTICAL  FORCE  ON  TOP  HALF  SEGMENT  OF  CABLE  DUE  TO  CURRENT 

C FWYIN)  ■ WEIGHT  OF  CABLE  SEGMENT  IN  AIR 

C FXIN)  • HORIZONTAL  FORCE  ON  CABLE  SEGMENT  DUE  TO  CURRENT 

C FXBUOV  • HORIZONTAL  FORCE  ON  BUOY  DUE  TO  CURRENT 

C FXCAB  - HORIZONTAL  FORCE  ON  TOP  HALF  SEGMENT  OF  CABLE  DUE  TO  CURRENT 
C FXMAST  • HORIZONTAL  FORCE  ON  ANTENNA  DUE  TO  WINO 

C FXPRBL  • HORIZONTAL  FORCE  ON  PARABOLIC  WEIGHTS  DUE  TO  CURRENT 

C FXTRPO  • HORIZONTAL  FORCE  ON  TRIPOD  DUE  TO  CURRENT 

C FYIN)  • TOTAL  VERTICAL  FORCE  ON  CABLE  SEGMENT 
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H • HEIGHT  OF  BUOY  ABOVE  WATERLINE 

KE(N)  - INVERSE  ELASTIC  MODULUS  FOR  CABLE  SEGMENT 

LOOPA  • ACCURACY  LOOP 
LOOPE  - ERROR  LOOP 

LORIG(N)  - ORIGINAL  UNSTRESSED  LENGTH  OF  CABLE  SEGMENT 
LT(N)  - LENGTH  OF  CABLE  SEGMENT  UNDER  TENSION 
LTRPD  - LENGTH  OF  TRIPOO  LEGS 

N ■ CABLE  SEGMENT  INDEX 

NL(N)  • NUMBER  OF  PARALLEL  LINKS  IN  A CABLE  SEGMENT 
NS  - NUMBER  OF  CABLE  SEGMENTS 

RHOAIR  - DENSITY  OF  AIR 
RHOWAT  - DENSITY  OF  WATER 

RXIN)  • HORUONTAL  RESULTANT  FORCE  ON  A CABLE  SEGMENT 
RYCNI  • VERTICAL  RESULTANT  FORCE  ON  A CABLE  SEGMENT 

SINTH(N)  - SINE  OF  ANGLE  BETWEEN  CABLE  SEGMENT  AND  VERTICAL  AT  THE 
N-TH  NODE 

T(N)  - TENSION  IN  CABLE  SEGMENT 

TBUOY  - THICKNESS  OF  BUOY 

TPRBL  - THICKNESS  OF  PARABOLIC  WEIGHTS 

VAIR  • WIND  VELOCITY 
VSBUOY  ■ SUBMERGED  VOLUME  OF  BUOY 
VTBUOY  • TOTAL  VOLUME  OF  BUOY 
VWAT  - VELOCITY  OF  WATER  CURRENT 

WCAB(N)  • WEIGHT  PER  FOOT  OF  CABLE  SEGMENT  IN  AIR 

WE(N)  - WEIGHT  OF  ELEMENT  AT  NODE  (N>  IN  WATER 

WTSYST  • TOTAL  WEIGHT  OF  BUOY  SYSTEM  W/O  CABLE  AND  ANCHOR 

X<N>  • HORIZONTAL  COORDINATE  OF  NODE  (N> 

XA  ■ XORIZONTAL  COORDINATE  OF  CABLE  END  AT  ANCHOR 

Y(N»  - VERTICAL  COORDINATE  OF  NODE  IN) 

YA  • VERTICAL  COORDINATE  OF  CABLE  END  AT  ANCHOR 
YCAB  - DESIRED  VERTICAL  COORDINATE  OF  TOP  END  OF  CABLE 

REAL  LOR IGI  20 1 *LT ( 20 ) #KE  I 20 ) *NL < 20 ) *LTRPD»MUWAT 
DIMENSION  DI A C 20 ) *D I 20 ) tWCAB ( 20 ) tCDCAB ( 20 > tCONST ( 20 ) • 

IDE  120) »WE(20) «T (20) .RX(20) tRY(20 ) #X( 20) »Y(20 ) .FX(20) »FY( 20) .X0(20 ) 
2 • YO ( 20 ) #FBY ( 20 ) *FCY ( 20 ) #FWY ( 20 ) t SI NTH ( 20 ) *COSTH ( 20 ) 

DATA  SINTH/20»0«0/*COSTH/20*1«0/ 

DATA  PI «DELTA/3* 14 159* 1000*0/ 

DATA  XA/0.0/ 

DATA  DELE »DELXY/0«0001 *0*1/ 

READ  IN  DENSITY  OF  WATER  AND  AIR  IN  LB.SEC««*2/FT*»*A*  WTSYST  IN 
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IUm  va  in  ft* 

9 I READ 11*101) RHOWAT  *RH0A I R • WTSYST • YA • NS 

101  FORMAT  UP10*S*1 9) 

READ  IN  LORIG  IN  FEET*  OIAMETER  IN  INCHES*  WCAB  IN  POUNDS*  WE  IN 
POUNDS 

DO  IS  N«1*NS 

IS  READI1*102)L0AI0IN)*DIAINI*WCAB(N)*C0CAB(NI*C0NST(NI*KE(NI*NL(NI* 

1DEIN)*WE!N) 

102  FORMAT! SF1 0*9 *F9*0*F 10*9 *F9 *2) 

READ (1*109 1 CDBUOY *CDPRBL *CDTRP0 

109  FORMAT (9F9.9) 

WRITE ( 3*901 ) 

WRITE! 9*902) 

902  FORMAT (Til** PROPERTIES  OF  THE  CABLE*///) 

DO  10  N«1 *NS 
WRITE! 9*909)N 

309  FORMAT! Til •• CABLE  SEGMENT* *12/) 

10  WRITE! 9*90*)LOAIG(N)*OIAIM)*MCAGIN)*CDCAB(N)*CONST!N)»KE(N)*NL(N) 
90S  FORMAT !T1S** ORIGINAL  LENGTH* *F10*9**  FEET* /T1S* ‘ORIGINAL  DIAMETER* 
1 *F10*9*  • INCHES*/TU**WEIGHT'*F10*9**  POUNDS  PER  FOOT */T IS* ‘DRAG  C 
20EFFICIENT**F10*9/T16* ‘ELONGATION  CONSTANT* *F10*9/T1S»* INVERSE  ELA 
3STIC  MODULUS' *F10*9/TlS» 'NUMBER  OF  PARALLEL  LINKS* *F4*0// ) 

J-NS-l 

DO  17  N*1*NS 
LT(N)»LORtG!N) 

DIA(N)«0IA!N)/12* 

17  D(N)«DIA(N1 

INITIALIZE  SYSTEM  DIMENSIONS 

TBUOY-32./12* 

DBUOY-96./12, 

DTRPO« 1*91/12* 

LTRPO-9S./12* 

TPRBL-4.4S7S/12. 

DPRBL-1. 

WRITE (9*901 ) 

301  FORMAT (1H1) 

WRITE (9*912) 

912  FORMAT (Til* • PROPERT I ES  OF  THE  BUOY*///) 

WRITEI9*919)TBUOY*DBUOY*TPRBL*DPRBL*DTRPD*LTRPD*WTSYST 
313  FORMAT  IT IS •• BUOY  THICKNESS* *F10*9* • FEET */TiS* ‘BUOY  DIAMETER* *F10. 
19*'  FEET */TlS*' PARABOLIC  WEIGHT  THICKNESS* *F10*9* ' FEET* /TiS* • PARA 
2BOLIC  WEIGHT  DIAMETER* *F10*9* • FEET* /TIS* 'DIAMETER  OF  TRIPOD  LEGS* 
9*F10*9* ' FEET* /TIS* 'LENGTH  OF  TRIPOO  LEGS* *P10*9* * FEET* /TIS* *WE IG 
AMT  OF  CABLELESS  BUOY  SYSTEM* *F10«4* • POUNDS'# IX) 

COMPUTE  TOTAL  VOLUME  AND  CROSS-SECTIONAL  AREA  OF  BUOY 

VTBU0Y»3*PPI<*TBU0Y**2 
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ATBUOY"2»*TBUOY*DBUOY/3« 

A-ARSIN07./98.) 

B«0«S2U988 

READ  IN  DEPTH  IN  FEET*  VELOCITY  OF  CURRENT  ANO  WIND  IN  KNOTS 

104  READ! 1 * ICS \DEPTH*VWAT *VA1R 
103  FORMAT! 3F10«2) 

NEGATIVE  NUMBER  FOR  DEPTH  ALLOWS  FOR  COMPLETE  RECYCLE  OF  INPUTS* 
BLANK  CARD  TERMINATES  PROGRAM 

IF (DEPTH) 31 *100*21 
21  VWAT« 1 • 66694* VWAT 
VAIR«1«68694*VAIR 
VWATSQ«VWAT**2 
VSBUOY-VTBUOY/2. 

ASBUOY-ATBUOY/2. 

YCAB-DEPTH  - 9.3 

LOOPE-O 

LOOPA-O 

MUWAT-RHOWAT  * VWATSO  / 4* 

COMPUTE  CONSTANT  FORCES  (POSITION  INDEPENDENT) 

FXMAST*0«012911#ABS( VAIR ) *VAIR 

FXTRPD-0*  3*CDTRPD«RH0WAT*DTRPD*LTRPD*  ( COS  ( A ) *«S*2*  ( COS  ( B ) *«3<*S  I N ( B 
1 I **3*C0S ( A ) **3 I ) *VWATSQ 

FXPRBLa2**0#3*RHOWAT«CDPRBL*2**TPRBL#DPRBL*VWATSQ/3. 

FBTRPD-3#*RH0WAT*32*174«PI*DTRPD*»2*LTRPD/4» 

FBPRBL*2 • *RHOWAT  *32 • 174*P 1 *TPRBL*DPRBL**2/16 • 

FBCAB*RHOWAT*32 • 174*P I *0 I A ( NS ) **2*L0R I G ( NS ) *NL ( NS ) /4 . 

DO  13  N«1*J 

FWY(N ) ■ (LORIG  ( N ) *WCAB < N ) *NL  (N  l+LOR IG( N+l  )*WCAB(N+1  )*NL(N*'l ) I / 2«  + 
1WE(N) 

13  FBY (N )aRHOWAT*32»174#PI*(DIA(N ) **2*LORIG(N)*NL (N ) ♦DIACN+l  1**2* 
1L0RIG(N+1  )#NL(N+1)  ) / 0 • 

COMPUTE  POSITION  DEPENDENT  FORCES 

FCYCAB«MUWAT*CDCAB ( NS  > «D ( NS  > »LT ( NS ) «S INTH ( NS ) *COSTH (NS ) **2*NL ( NS  > 
18  LEAP*1 

FXBU0Y»0«3*RH0WAT*CDBU0Y*ASBU0Y*VWATS0 
FXCABaMUWAT*CDCAB ( NS ) *D ( NS ) *LT ( NS ) *C0STH ( NS ) **3*NL ( NS ) 

FX  ( NS ) •FXMAST^FXBUOY'fFXTRPD+FXPRBL+FXCAB 
FBBUOY*RHOWA  T*32 • 1 74*VSBUOY 

FY  ( NS ) -FBBUOY-fFBPRBL-f  FBTRPD+FBCAB-WTSYST-FCYCAB 
DO  11  Nal • J 

FX  (N  ) aMUWAT*(CDCAB(  N)*D(N)  *LT  ( N ) *C0STH( N )**3*NL (N ) ♦(COCAS  < N+1 ) #D  (N 
1*-1 ) *LT  ( N<*1  )*NL ( N*-l ) +DE  (N  ) ) #COSTH ( N+l ) **3 ) 

FCY ( N ) aMUWAT* ( CDCAB ( N ) *D ( N ) *LT ( N ) *SI NTH ( N ) *C0STH ( N ) **2*NL ( N ) ♦CDCAB 
i(N‘*l)*D(N4l)*LT(N+l)«SINTH(N-*l)#C0STH(N-*l)**2*NL(N+l)  > 

11  FY(N>aFBY(N)-FWY(N)-FCY<N> 
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COMPUTE  RESULTANT  FORCES 

9 RXtNSI-FXCNS) 

RY(NS>>PY(NS) 

00  9 N«1*J 
N«JH*»1 

RX(N)*FX(N)  ♦ RX(N+1) 

9 RY(N)»PV(N)  ♦ RYtN+l) 

COMPUTE  SEGMENT  TENSIONS  AND  LENGTHS  UNDER  TENSION 
DO  7 N-l.NS 

T C N I •SORT ( RX ( N ) »*2*RY ( N ) #*2 1 

LT (N)*LORIG(N )*( 1*0*C0NST(N )+T (N)*KE (N) /NL(N ) ) 

IF(LTIN) /L0R1G(N)*LT*  1.1)  GO  TO  • 

D(N)»DtA(N)*SQRT(LORIG(NI/LT(N)l 
8 SINTH(N)*RX(N)/T(N) 

C0STH(NI«RY(N)/T(N> 

COMPUTE  NODE  COORDINATES 

GO  TO  (4t5)»N 

S X ( N ) *L T ( N ) *S I NTH ( N ) +X ( N*1 > 

V I N I »LT ( N ) *C0STH ( N I ♦Y ( N«1 I 
GO  TO  7 

4 X ( N ) «LT ( N ) *S I NTH ( N I ♦XA 
YIN)»LT(N)«C0STH(NI4>YA 
7 CONTINUE 
L00PE*L00PE^1 

COMPUTE  AND  COMPARE  ERROR  FUNCTION 
E*tYCAB-Y<NS)>««2 

IF(E*GT*OELE)  GO  TO  < 1Y»50* 15 I »LEAP 
DO  5S  N-l.J 

IF! ABS(X(N)*X0(N) I •GT*DELXY*OR«ABSl Y (N)»Y0I N ) ) »GT *DELXY)  GO  TO  57 
55  CONTINUE 

COORDINATE  ACCURACY  SATISFIED  - WRITE  RESULTS 

WRITE < JtSOl  I 
WRITEOtSOSI 

505  FORMAT! Til ••GUOY  AND  CABLE  CONFIGURATION  UNDER  CURRENT  AND  WlND'/T 
195* • CONDITIONS*///) 

WRITE (9»908>Xt NS >»DEPTH»VWAT*RH0WAT»VAIR»RH0A1R 
908  FORMAT! TIG •• BUOY  HORIZONTAL  DISPLACEMENT  IS*»F10.9»*  FEET'/TIG* 'OE 
1PTH  OF  WATER  IS**F10«9>*  FEET' /T 16#' VELOCITY  OF  WATER  IS*»F10.5«* 
ZFEET  PER  SECOND* /T 18* • DENSITY  OF  WATER  IS'tFlO.5**  SLUGS  PER  CU«FT 
9* */Tl8* 'VELOCITY  OF  AIR  IS*»F10«5»*  FEET  PER  SEC0ND*/T16»*DENSITY 
80F  AIR  IS* •F10»9«*  SLUGS  PER  CU.FT.*/) 

H-12,*H 

WR I TE  < 9 • 914 1 FXMAST .FXBUOY • FXTRPD  tFXPRBL • H 
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914  FORMAT I T16  » ' ANTENNA  DRAG  IS*.F10.9«*  POUNDS* /Tit ** BUOY  DRAG  IS*. FI 

10.9. '  POUNOS*/TU.* TRIPOD  DRAG  IS'.F10«9.'  POUNDS '/T 16. 'PARABOLIC 
2 WEIGHT  DRAG  IS'.FIO.S.'  POUNDS' /TIG. 'WATERLINE  IS'.F9#2.»  INCHES  F 
IRON  TOP  OF  BUOY*///) 

DO  S*  N*i .NS 
TIN)  •TIN) /NUN) 

99  WRITEI9»911)N»XtN).YIN).N.TlN)»N.FXIN) 

911  FORMAT  I TIG. 'COORDINATES  FOR  NODE  I*.  1 2.')  ARE  X«'.F10«9.'  FEET. 
1Y«* .F10.9.'  FEET'/T21. 'TENSION  IN  SEGMENT  l'.I2.')  IS'.FIO.S.'  POU 
2NDS'/T21. 'HORIZONTAL  FORCE  AT  NODE  (*.I2.«)  DUE  TO  CURRENT  IS'.FIO 

9.9. '  POUNDS'//) 

WRITE  1 9.909) tOOPE.LOOPA.E .DELTA 

909  FORMAT! Til. 'NUMBER  OF  ERROR  LOOPS  •' .19. T49* 'NUMBER  OF  ITERATION  L 
IOOPS  •'.19 //Til. 'ERROR  FUNCTION  •' .E14.7/T11 . 'DELTA  ■•.E14.7) 

GO  TO  104 

COORDINATES  NOT  WITHIN  DESIRED  ACCURACY  - RECORD  NODE  COORDINATES 
FOR  COMPARISON 

97  DO  97  N*1.J 
X0IN)*XIN) 

97  Y0IN)«YIN) 

RECOMPUTE  SUBMERGED  BUOY  VOLUME  AND  CROSS-SECTIONAL  AREA 

FCYCAB«MUWAT*CDCAB  I NS ) «D  I NS ) *LT  I NS  > "S I NTH  I NS  > "COSTH  INS)  «*2*NL  ( NS  ) 
FBBU0Y«FY INS ) +WTSYST+FCYCAB-FBPRBL-FBTRPD-FBCAB 
VSBUOY»FBBUOY/IRHOWAT«92.174) 

H-SORTI IVTBUOY-VSBUOY) /I4.»P1 ) ) 

ASBUOY • ATBUOY-8 • •HwSOR  T ( 9 • «H ) / 9 • 

DELTA-FYINS) 

LOOP A* LOOP A ♦ 1 

RECOMPUTE  POSITION  DEPENDENT  FORCES 
GO  TO  18 

CHECK  IF  ERROR  FUNCTION  HAS  DECREASED 
90  IFIE.LT.EP)  GO  TO  20 

IF  ERROR  FUNCTION  HAS  NOT  DECREASED.  DECREASE  CONVERGENCE  FACTOR 

OELTA-DELTA/2. 

12  DELT-OELTA/SORTIEP) 

CHANGE  VERTICAL  FORCE  TO  MOVE  END  COORDINATE  TO  DESIRED  VERTICAL 
POSITION 

FY I NS ) *FYP* I YCAB-YP ) *DELT 
IF  I FY I NS) .NE.FYP ) GO  TO  3 

NO  CHANGE  IN  FORCES  - TERMINATE  PROGRAM 
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LEAP-) 

GO  TO  ) 

IS  WRlTEOtSOl) 
tfRITE(9*S07) 

307  FORMAT (Til* 'PROBLEM  IS  NOT  COMPLETE*  DELTA  HAS  BECOME  TOO  SMALL  TO 
1 CHANGE  THE  IMAGINARY  REACTIONS*) 

GO  TO  100 

19  GO  TO  (22*2)>  *LOOPE 

SAVE  NODE  COORDINATES  FOR  COMPARISON 

22  DO  14  M-l*j 
XOIID-XIN) 

14  YO(N)-Y(N) 

2)  LEAP-2 

[C  SAVE  ERROR  FUNCTION*  CABLE  END  COORDINATE  ANO  FORCE 

| 20  EP-E 

YP-YCNS) 

[ FYP-FY(NS) 

k GO  TO  12 

I 100  CALL  EXIT 
ENO 


